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Abstract 



Generalized action invariants are identified for various models of 



^^ I drift wave turbulence in the presence of the mean shear flow. It is 

o 

, shown that the wave kinetic equation describing the interaction of 

the small scale turbulence and large scale shear flow can be naturally 
, writen in terms of these invariants. Unlike the wave energy, which is 

jj^ conserved as a sum of small- and large- scale components, the gener- 

alized action invariant is shown to correspond to a quantity which is 
conserved for the small scale component alone. This invariant can be 
used to construct canonical variables leading to a different definition 
of the wave action ( as compared to the case without shear flow). It 
is suggested that these new canonical action variables form a natu- 



ral basis for the description of the drift wave turbulence with a mean 
shear flow. 

The dynamics of the small scale turbulence in the presence of a mean 
shear flow is a problem of a great interest for plasmas and geostrophic fluids. 
It is believed that the nonlinear energy transfer from small to large length 
scale component (inverse cascade [^ ) is a cause of a spontaneous generation 
and sustainment of coherent large structures, e.g. zonal flows in atmospheres, 
ocean and plasmas 0. In the few past years it has been suggested p, ^, ^ 
H) 0) H) B 113 that the large scale flow band structures (zonal flows) play 
an important role in regulating and suppressing the anomalous transport in 
magnetic conflnement systems. 

In the simplest form, the generation of plasma flow by turbulence can be 
described by the energy conservation relation (Poynting theorem) averaged 
over small scale fluctuations |^ . A generalization of this approach is a WKB 
type wave kinetic equation for the quanta density of small scale fluctuations 
that is conserved along the rays. This method was originally proposed in 
Ref. 11 to describe the interaction of high frequency plasmons (Langmuir 
waves) with low frequency ion sound perturbations. It is widely used also in 
general fluid dynamics fir 



In studies of drift wave dynamics, it has been naturally assumed ||13| , |I4 
that the relevant quantity that is conserved in the presence of slow variations 
is the drift-wave action density. It is well known |]T5|, that the standard wave 



action variables Ck associated with the number of wave quanta Uk, rik = 
\Ck\ = Ek/ujk, where Ek is the wave energy, and Uk is the wave frequency, 
is a basis for Hamiltonian form of the wave-wave interaction equations. It 
has been noted in Refs. 16,17 that the normal variables used to describe 
self-interaction between small scale fluctuations without the shear flow are 
modified by the flow and may not be suitable for a system with a mean 
flow. Thus, in the presence of a shear flow a new form of canonical variables 
and associated action invariant have to be identified. On other hand, it has 
been pointedfl^ that the conserved action-like quantity (pseudo-action) is 



different from the usual definition of the wave action defined as the ratio 
of the wave energy to the wave frequency. The latter definition is also fails 
when there are no oscillating eigenmodes such as in ideal fluid, so that an 
alternative definition of the action-like integral is required ||19||. 

It is important to realize that the natural form of the three-wave in- 
teraction equations for the drift-waves does not have Hamiltonian structure 
[^. These equations can be transformed, however, to a Hamiltonian form 
via an asymptotic variable transformation. Such a tranformation yielding a 
Hamiltonian form for the drift and Rossby waves has been found in Refs. 
20,21. There are several possible forms for such a transformation. In Refs. 
17,20,21 it is based on the conserved energy integral that leads to the stan- 
dard definition of the wave action. For drift-wave+zonal flow systems small 
scales are modulated by larger scale shear flows so that energy in the small 



scale component is not conserved. Thus, the canonical Hamiltonian variables 
constructed from energy conservation are not suitable for description of the 
drift waves in the presence of a mean flow. 

In this paper, we derive the WKB type wave kinetic equation that de- 
scribes the conservation (along the rays) of an action like invariant of the drift 
wave turbulence with slowly varying parameters due to the mean sheared 
flow. We demonstrate that the relevant action-like integral corresponds to 
the quantity conserved for the small scale component alone. We show that 
the structure of the action integral is determined by the structure of the 
matrix element describing the interaction of the small scale and large scale 
component. We discuss how the canonical variables corresponding to such a 
pseudo-action invariant can be constructed. 

The scale separation between the small scale turbulence and the large 
scale motions is an essential property of drift-wave-|-zonal flow systems that 
is commonly used (jT^], |T^, [T^, |T^, |I^, ^, ^ to simplify the analysis. Though, 



the scale separation is often observed experimentally and in computer sim- 



ulation, it may be less pronounced in other cases[0]. In our present paper, 
we substantially rely on the multiscale expansion, so our results are valid, 
strictly speaking, only in the case when there is such a scale separation. More 
general approach avoiding the scale separation assumption, namely the renor- 
mahzation group, is possible p5[) but it is beyond the scope of the present 



paper. 



We consider a generic case of the drift wave equation in the form 

-^ + iuJk(t)k + / d^pLp,k^p(t)p(t)k-p = 0, (1) 

where Uk = u>{k) is the frequency of the hnear mode with a wavector k, and 
may include an imaginary part corresponding to the wave grow and decay. 

In the spirit of the scale separation we represent the field into the large- 
scale (p^ and small-scale (j)^ components; 0^ = outside a shell | k |< e ^ 1, 
0^ = for I k |< e. 

Assuming that the self-interaction of small-scale fields is small compared 
to the interaction with the mean flow||T^ we write from ([^) the following 



equation for the small-scale fluctuations 

+ iujk^> + I (fpLp^k-p<P^(p^-p = 0. (2) 



^r 



dt 

To derive the equation for the evolution of the wave spectrum we multiply 
equation (0) by (j)^, and then add it with a similar equation obtained by 
reversing k and k , yielding 

Ol (<^fc </'fc')+* ^^^ + '^k') 0fc 0fc'+0fc' J d'^vLp,k-p(t)p(t)l-p+4>l J d^pLp^k'-p<Pp(Pt'^p = 0- 

(3) 
The small-scale turbulence is described by the spectral function (Wigner 
function) /^(x, t), defined as follows 

J d\ {4>\+^4>l) exp(«q ■ x) = 4(x, t). (4) 



The slow time and spatial dependence in Ik{'x,t) corresponds to modu- 
lations with a "slow" wavevector, q ^ k. Angle brackets in (^) stand for 
ensemble average, which is equivalent to a time average with appropriate 
ergodic assumptions. 

The equation for /^(x, t) is derived from (^) by averaging it over fast 
scales and by taking the Fourier transform over the slow variable x. Setting 
k = — k + q and applying the operator / c/^gexp(iq • x) we obtain 

—h(x,t)+i J (fqexp(iq-x.)(iJk + uj^k+q){(l)k4>^k+q) + 5*1 + 5*2 = 0, (5) 

•^i = y y d^pd'^qexp{iq ■ x) (^(p^k+g^Pk-p) Lp^k-p<Pp, (6) 

S'i = j j d^pd^qex^iiq, ■ x) {(Zk+q-p4>l) Lp,-k+q~p(pp ■ (7) 

The second term in (|^) gives 

i j d^q exp(iq ■ x) {ujk + uj^k+g) {(Pk'P-k+g) = "^ " g^^''^'^^ ^) ~ 27fc4, (8) 

where 7^ is the linear growth rate, and only the real part of the frequency is 
presumed for Uk on the right hand side of this equation. 

The ensemble average in 5*1 can be transformed by using the inverse of 



'^-k+g4>k-p) — \4>k-p4>-(k-p)+q- 



p) = J dV/fc_p(a;')exp(-i(q-p)-x'). (9) 



By using (^ and expanding in p ^ k the expression for Si is transformed 
to 

Si = J rfVexp(zp ■ x)Lp^k-p ( 4(x) - p • ^^ J (pp . (10) 

Similarly, by using the identity analogous to (^ and expanding the interac- 
tion coefficient Lp^k_p in p <^ k, we transform 5*2 to the form 

S2 = J J d^p(fq exp(iq • x) ( Lp^_k + (q - p) ■ -qT^ ) '^p 
X / ci^x exp(— ^(q — p)-x )/fe(x ) 

d'^pexp{ip ■ x.)Lp^^k(j)< -i J d^pexp{ip ■ x) ^^'^^ ■ -^0|ll) 

Equations (^-|Tl|) define a particular form of the transport equation for 
/fc(x, t) for a given interaction coefficient -Z^^^,'. 

In this paper, we consider two different models for drift waves in a magne- 
tized plasma: the standard Hasegawa-Mima equation and a slab-like model 
for drift waves in a sheared magnetic field. The latter is similar to the stan- 
dard Hasegawa-Mima equation with a modified plasma response to the slow 
modulations of the electrostatic potential. Such slow modes correspond to 
^11 — s> 0, so that the slow part of the potential does not follow Boltzmann 
distribution. [Note that zonal flows |jlO| {m = n = 0) are such slow modes 
with k\\ = 0.] As a result, the convective term appears in the lowest order, 
contrary to the case of the Hasegawa-Mima equation where such term is due 
to the polarization drift. Appropriate equation for the drift wave dynamics 



in presence of a mean flow (neglecting the self-interaction) has the form []T3 



(|,v..v)£f.v..v^-4|.v„.v)vi^^o. a, 

where Vq = cbxV0/-Bo is the mean flow velocity. This equation can be 
written in the form (||) with cjfc = k • V*/(l + fc^p^) ^^^ 



_ c b ■ ki X k2 (^ ,22 

'''''' ~ ~s;i + (ki + k2)V 



1 + fc^V' ) . (13) 



^s 



From (|5|-p!ID and ([131) we obtain 



^4(x. ,) + _ (., H- k . V„) . - - - ( (^^^j ■ 5i^4(H- *V^)'^ = 0. 

(14) 

This equation can be written in the form of a conservation law for the in- 
variant Nk = 4(1 + k^p'^y, 

5^7V.(x.«, + -(.. + k.V„).^--(k.Vo)-A'.^0. (16) 



By direct evaluation from (|T2D, it can be easily shown that the quantity 
N = jd'k (02 + 2pl{V^4>f + p'siyW) , (16) 

corresponding to A^^ in (17), is conserved as an integral over the small-scale 
part of the spectrum. In (p!6D is the normalized potential of the small scale 
fluctuations. This property distinguishes Nk from any other combination of 
the energy and enstrophy which are conserved only as a sum of contributions 
from the small and long scale parts of the spectrum[^ . 



A different expression for the action-like invariant is obtained for the 
standard Hasegawa-Mima (H.M.) model with a mean flow 



m-^y4]....Vf-Aiy...)V' 



0. 



The appropriate interaction coefficient is 



c 2 b ■ ki X k2 

Ps 



\^2 ^l) ■ 



(17) 



(18) 



2i?o"n + (ki + k2)2p2 

In this case, from (5-11) and (|TB|) the transport equation for /^ takes the 
form 



m'' + 9k h + TTWs '' ■ ^ ' 



d 



k V 







d 



k'pl{l + k^pl)h 



0. 



(19) 



9x V(l + A;2ps)V ^k 
Obviously, this equation can be written in the form of the conservation law 



for the invariant N^ = hk^Ptil + A;Vs)0 El 
9 ^r 9 f k-Vo ,2 2\ ^ .. d 



^■^° -k'pl]-^N, = 0. 



k'pl) 



d 



(20) 



Similarly to the previous case, the invariant Nk corresponds to the integral 
of (p!7D conserved for the small scale component alone 



N 



d'k pIUv r4>f + pI{v] 



(21) 



Note that both invariants (|T^) and (^) are different from standard definition 
of the wave action |T3|, |1^. The difference between two forms of the action- 
like invariant (Eq. (|1^) and (0)) is due to a different form of the coupling 



matrix (Eq. ( |I3| ) and Eg. ([I8|) ) describing interaction of the small and large 
scale components. 

The procedure that we have described above can also be used to derive 
the action-like invariant for the two-dimensional motion of an incompressible 
fluid. In the latter case, there are no oscillating modes so that the standard 
definition of the action as a ratio of the wave energy to wave frequency is not 
applicable. The 2-D Euler equation has a form 

avl0 + Vo ■ VVl0 = 0, (22) 

where Vq is the velocity due to the mean flow. This equation can be written 
in the form (|^) with uJk = and the interaction coefficient 



_ b ■ ki X ka J 2 



kl (23) 



Using equations (5-11) and (ESf) we obtain the wave kinetic equation 



-N.i.,t)^-ik.V„).^--i,..V„).-N.^O, (24) 



where the wave-action A^^ = k h ||T9|. 



We summarize generalized wave action integrals for different models in 
the Table I. Note that the standard expression for the drift wave action 



defined as the ratio of the wave energy to the wave frequency is [O, 114 



nk = \ak\ = \(Pk\ = — , (25) 

where cu* = /c^K- Expression (^) should be compared with the first two 
lines in the Table. It is interesting to note that generalized action invariant 

10 



given by Eq. (|T^) coincides with the standard definition of the wave action 
( p5D for the poloidally symmetric case when the poloidal wave vector kg is 
constant {kg = const). 

Next we consider the self-interaction between small scales in the presence 
of the shear fiow and outline how the pseudo-action invariants can be used 
to construct the canonical variables for the latter case. For illustration, we 
consider the case of Hasegawa-Mima equation ([T7|) . We restore the self- 
interaction term given by Wk^k^M 

'' + iujk(j)k = / d'^ki(fk2Wk,kuk2S{k -ki- k2)(j)ki(j)k2, (26) 



dt 

This natural form of the three-wave interaction does not have standard 
Hamiltonian structure. This is refiected in the interaction coefficients Wk^kiM 
which do not have the required symmetry properties |T^. The only symme- 
tries in WkMM are of the type Vr*^.^^.^^ = Wk^M = W-kMM- Transfor- 
mation of ( P^D to normal canonical variables a^ was given in Refs. 19,20 (see 
also Ref. 16). It has the form |]17 



Qfe = gk(t>k + / (fki(fk2Gk,kuk2^{k -ki- k2)(t)kAk2- (28) 

In new variables the interaction coefficients Vk^kiM ^^^ 

Vk,ki,k2 = n \\9k\ Wk,ki,k2 + Ifi'fcil Wkl^k,k2 + Ifl'tel Wk2.ki,k) (29) 

^9k^9k2 9k ^ ' ' 



11 



These interaction coefficients Vk^kiM ^^^ have all symmetries required for 
Hamiltonian systems. The function gk can be chosen in a variety of ways. The 



standard approach [0, ^, |21| is to chose gk so that the energy in canonical 
variables takes the form E = ^ (PkukCikO'-k- Comparing it with the energy 
integral E = J £k {^^ + p^iV ±$y)) ,we find ^ gk = (l + plkl)/{kyf'\ 
This gives a standard expression for the wave action (p5[) . 

As discussed above, for the drift waves-zonal flow system the energy in 
the small scale component is not conserved, bur rather the total energy of 
drift waves + large scale zonal flows is constant. For this reason, the energy 
integral of the small scale component can not be used for introduction of 
canonical variables for self interaction of the small scale fluctuations. Con- 
trary to the energy, the integrals A^^ are conserved for small scale component. 
Choosing the function gk such as that the invariants (|T6|) or (^) are in the 
form Nk = ^ (Pkaka^ki^^ obtain Nk as canonical variables for drift waves in 
the presence of the mean shear flow. This automatically means that these 
invariants have a meaning of the generalized wave action invariant. Then, 
to account for the self-interactions in the presence of the background shear 
flow, the wave kinetic equation (Eq. (|15D or (^OD) should be modified with a 
source term Jk in the standard form||15|| 



Jk = All j (fki(fk2 X \VkMM? (Nk.Nk, - NkNk, - NkNk) 6{k - ki - /ca). 

(30) 
We have formulated a wave kinetic equation and determined a structure of 



12 



an appropriate adiabatic invariant for small scale turbulence in the presence 
of a mean flow. We have shown that the form of the matrix coefficient for 
the nonlocal coupling of the small scale fluctuations to the mean flow is 
crucial for the form of the adiabatic invariant. We have obtained adiabatic 
invariant A^^ = ^kk'^pli^ + k'^p'i) for the drift wave turbulence described by 
the Hasegawa-Mima equation and isomorphic Charney-Obukhov equation for 
Rossby waves; and the invariant Nk = /fc(l + k'^plY ^^ ^^e drift wave type 
turbulence in tokamaks such as TITG driven modes. [Note that the latter 



invariant reduces to the standard form [|13], Q for kg = const.] The pseudo- 
action invariants appear in the wave kinetic equation and correspond to the 
quantities that are conserved as integrals over the small scale part of the 
spectrum alone. This speciflc conservation property makes them suitable as 
canonical Hamiltonian variables for small scale turbulence in the presence of 
the shear flow. The wave action invariants and the kinetic equation derived 
here can be used to investigate nonlinear dynamics of drift waves and zonal 
flow in a tokamak. The method used in our work can be applied to derive 
generalized invariants for other models including the Rossby type waves in 
geostrophic fluids [|l^. 
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Table I: Generalized action invariants for different models 

Model Expression for the wave action 
Drift waves in a sheared field, Eq.(|12|) 4(1 + k'^pi)'^ 

Standard drift wave model, Eq.(|l3) Ikk'^p^l + k'^pl) 

2D Euler equation, Eq. (|2|) hk'^ 
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